Simple and efficient generation of gap solitons in Bose-Einstein condensates 
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We suggest an efficient method for generating matter-wave gap soiitons in a repulsive Bose- 
Einstein condensate, when the gap soliton is formed from a condensate cloud in a harmonic trap 
after turning on a one- dimensional optical lattice. We demonstrate numerically that this approach 
does not require preparing the initial atomic wavepacket in a specific state corresponding to the edge 
of the Brillouin zone of the spectrum, and losses that occur during the soliton generation process 
can be suppressed by an appropriate adiabatic switching of the optical lattice. 
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I. INTRODUCTION 

Recent experimental developments in the fields of non- 
linear optics and Bose-Einstein condensation provide 
unique opportunities for the studies of many nonlinear 
phenomena. Moreover, due to a formal similarity of 
the Gross-Pitaevskii equation for the dynamics of matter 
waves and the nonlinear Schrodinger equation for light 
propagating in a nonlinear Kerr medium, it is possible 
to interchange and transfer many interesting ideas be- 
tween these two fields. One of them is the generation of 
nonspreading localized wavepackets usually called optical 
solitons or matter-wave solitons [J. 

Bright solitons, which are quite common objects in 
different problems of nonlinear optics, have been only 
recently demonstrated in two independent experiments 
with matter waves 0, Q- In this case, the wavepacket 
spreading due to dispersion is compensated by an at- 
tractive interaction between atoms. Both optical and 
matter-wave solitons can also be formed due to an inter- 
play between nonlinearity and periodicity ^ 1^ - 

In optics, the presence of a periodically varying refrac- 
tive index of a medium results in a bandgap structure 
of the transmission spectrum of light. This subsequently 
affects the propagation of optical beams (pulses), which 
is determined by the diffraction (dispersion) curves of the 
transmission bands and forbidden gaps of the wave spec- 
trum. It was demonstrated that periodic structures can 
modify dramatically the wave diffraction 7j, and there- 
fore they introduce a novel way for controlling self-action 
of light in nonlinear periodic media JJ. In particular, 
if the refractive index decreases with the light intensity 
due to a nonlinearity-induced response of the material, 
the beam normally experiences broadening due to self- 
defocusing. However, in periodic photonic structures, 
the same type of defocusing nonlinearity can lead to the 
beam localization 10] . For deep periodic modulations, 
this effect should result in the formation of discrete soli- 
tons. Indeed, as was shown both numerically and exper- 



imentally in Ref. , increasing the depth of the refrac- 
tive index modulation leads to a sharp crossover between 
the continuous and discrete properties of the lattice. It 
is manifested in a specific type of light localization and 
gap-soliton generation which become possible for a self- 
defocusing medium only in case of a discrete model. 

Analogous effect of the gap-soliton generation in the 
regime of repulsive interaction was observed for the Bose- 
Einstein condensates. This first experimental observa- 
tion ^2 was reported for *^Rb atoms in a weak optical 
lattice. In this case, the anomalous dispersion of the 
atomic cloud is realized at the edge of the Brillouin zone; 
if the coherent atomic wave packet is prepared in the 
corresponding state, a matter-wave gap soliton is formed 
inside the spectral gap. However, it took a consider- 
able experimental effort to generate the initial state of 
the Bose-Einstein condensate in a proper way, by intro- 
ducingan acceleration to the optical lattice and other 
trickslHlli. 

In this paper, we reveal that the matter-wave gap soli- 
tons can be generated by a simpler and straightforward 
method that allows to create stable, long-lived gap soli- 
tons of a repulsive Bose-Einstein condensate in a one- 
dimensional optical lattice. In our scenario, the gap soli- 
ton is formed from an initial ground state of a harmonic 
trap after turning on the optical lattice, due to the cloud 
relaxation rather than as a result of preparing the ini- 
tial atomic wavepacket at a certain point of the band 
gap structure. In addition, we demonstrate that losses 
during this process can be suppressed by an appropriate 
adiabatic switching of the optical lattice. 

The paper is organized as follows. In Sec. II we intro- 
duce our theoretical model. Then, in Sec. Ill we present 
two methods for generating the matter-wave gap soli- 
tons and demonstrate their efficiency in numerical simu- 
lations. The coupled-mode theory for describing the self- 
trapping mechanism is introduced in Sec. IV where we 
discuss the essential physics of the localization process. 
Finally, Sec. V concludes the paper. 
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II. THEORETICAL MODEL 

We analyze our matter-wave system in the framework 
of the three-dimensional mean-field Gross-Pitaevskii 
equation for the macroscopic wave function ^(r,t), writ- 
ten in the standard form, 

at \ 2m m J 

where for the repulsive condensate the scattering length 
is positive (a > 0) . We assume that the external potential 
U{r,t) consists of two parts, a harmonic trap and one- 
dimensional optical lattice, i.e. 

t/(r,t) = e{t)sin' (^) + | [co'ig^ + f{t)ulz^] , 

where d — A/2 is the period of the optical lattice poten- 
tial, p — [x^ + y'^Y^'^ is the two-dimensional radial co- 
ordinate, and the cylindrically-symmetric harmonic trap 
is characterized by the longitudinal and transverse fre- 
quencies, ojz and wj., respectively. We also introduce 
the depth e and period A of an optical lattice. Temporal 
variations of the potential are described by two functions, 
fi{t) and /2(i), which account for the adiabatic processes 
of turning off the longitudinal part of the harmonic trap 
and turning on the optical lattice over the stationary con- 
densate. More details about this tuning will follow in the 
next section. 

In our numerical calculations, we use both three- 
dimensional model and (more often) quasi-one- 
dimensional reduced model, in which we assume 
that the transverse harmonic confinement is so tight 
that no transverse excitations occur [13 . ITsI l . As a result 
of this assumption, the three-dimensional wave function 
resembles the transverse ground state of a harmonic 
potential in the transverse cross-section, and the non- 
linear dynamics can be described by a one-dimensional 
Gross-Pitaevskii equation after rescaling the nonlinear 
coefficient, aio — {mujj_/2'Kh)a. 

III. NUMERICAL RESULTS AND 
DISCUSSIONS 

In our numerical simulations, we suggest and analyze 
systematically two methods for generating matter-wave 
gap solitons in the repulsive Bose-Einstein condensates 
loaded in a shallow one-dimensional optical lattice. 

A. The first method 

In the first method, we start from the ^^Rb conden- 
sate (atomic mass 86.9 amu and the scattering length 
5.3 X 10~^ m) as an eigenstate of a cylindrically symmet- 
ric harmonic trap with the frequencies cozo = 27r x 40Hz 
and L0± = 2TT X SOHz. Then, we instantaneously turn off 




FIG. 1: (color online) Generation of a matter-wave gap soli- 
ton. Shown are the density snapshots for (a,b) = 200 atoms 
at the input {t — 0) and after expansion {t = 0.4s), respec- 
tively, and (c) for A'^ — 500 atoms, after expansion at t = 1.4s. 
The lattice parameters are: e — 2.3£'rccoii and period 15/im. 

the longitudinal trapping potential leaving the perpen- 
dicular harmonic potential on and simultaneously switch 
on, in the same direction, a one-dimensional optical lat- 
tice with the spatial period A = 30/im. For this lattice we 
define the recoil energy E'rocoii = ?i^(27r/A)^/(2m). Such 
a lattice can be created by interfering two laser beams 
propagating at a small angle |^ 0| . For the wavelength 
of 783 nm, used in Ref. [TjI, this angle would be equal 
to 2.3 degrees. After an abrupt change of the potential, 
a quasi-one-dimensional dynamics takes place, in which 
atoms slowly tunnel to the neighboring lattice sites, while 
the confinement in the direction perpendicular to the op- 
tical lattice remains unaffected. As is shown in Figs.^^a- 
c), depending on the initial number of atoms, we ob- 
serve asymptotically either expansion of the condensate 
over more and more lattice sites [see Fig. Wi^)] or con- 
densate stabilization accompanied by the formation of a 
spatially- localized state occupying just a few lattice sites 
[see Fig.^c)]. We checked that there exists a threshold 
in the initial number of atoms (with all the rest of the 
system parameters fixed) above which a stable spatially 
localized state is generated in the optical lattice. 

These results are summarized in Figs.|5Ja-c), where we 
show the profile of the soliton wavefunction and the cor- 
responding dynamics in the Fourier space. In the Fourier 
space, a clear localization at the edges of the first band 
occurs. To understand importance and specific features 
of this localization process, we notice that the evolution 
of a condensate in a one-dimensional optical lattice can 
be described within the effective mass approximation 

dtp d'^ip 47racff?i^, |2, 

'''m=-2^d^ + ^^^'^^^^ 

where ip(z,t) is an envelope of the condensate wavefunc- 
tion at the edge of the Brillouin zone of the spectrum. 
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FIG. 2: (color online) (a) Temporal evolution (from left to 
right) of the axial density of the condensate [see Fig.^c)]. (b) 
Evolution shown in the Fourier space, (c) Temporal variation 
of the depth of the optical lattice (cf. Fig. 5 below). 




FIG. 3: (color online) Wavefunction of the matter-wave gap 
soliton of Fig. Wic), as predicted by the one-dimensional 
Gross-Pitaevskii equation (dashed) and the full three- 
dimensional model (solid); Dotted line shows the lattice po- 
tential. 



The localized state can be supported only if the sign 
of the effective mass is opposite to the sign of the 
scattering length Ucs- Since for a repulsive condensate 
flcff > 0, the effective mass in this case has to be neg- 
ative (anomalous diffraction). This condition is indeed 
fulfilled for the matter waves close to the edge of the 
Brillouin zone. In Sec. V we present a simple model that 
allows to describe this kind of the localization process. 
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FIG. 4: (color online) Evolution of the peak amplitude of the 
atomic cloud in the lattice for (a) A'' = 200 atoms (solid), 
(b,c) A'^ = 500 atoms in the one-dimensional (dashed) and 
three-dimensional (dotted, rescaled) models, respectively. 



In Fig. 121 we show the wavefunction profile for the 
case of iV = 500 atoms a.t t = 1.4s, when the gap soliton 
is generated, as shown in Fig. ^c). The phase struc- 
ture of the resulting state corresponds to the edge of the 
first band of the Brillouin zone addition, we 

compare the results obtained from the integration of the 
one- and full three-dimensional Gross-Pitaevskii models, 
and notice that the condensate dynamics looks practi- 
cally the same in both the cases. With a good accuracy, 
we can therefore assume that the condensate remains in 
the ground state of the transverse harmonic potential 
that is left on all the time. Also, in Fig. 21 we show the 
evolution of the peak amplitude of the atomic cloud that 
fully supports our observation. As mentioned above, for 
N — 200 atoms we predict overall spreading of the wave 
packet and the peak amplitude vanishes gradually. For 
larger number of the input atoms (e.g., N = 500), the nu- 
merical simulations of both one- and three-dimensional 
models show the stabilization effect and both approaches 
predict almost the same final distribution of atoms on a 
few neighboring lattice sites. In earlier studies [T^ ITU, 
the generated gap soliton was unstable due to the excita- 
tion of transverse modes of the harmonic trap. Here, we 
suggest applying the optical lattice with the period com- 
parable to the transverse extent of the condensate. Such 
a configuration (i.e. a sparser lattice) has much more 
favorable ratio of the energies of the transverse and lon- 
gitudinal modes, and the suppression of the instability is 
observed. 



B. The second method 

In the loading scenario discussed above, the key idea 
was to replace suddenly the harmonic trap in the lon- 
gitudinal direction with a shallow periodic potential of 
a one-dimensional optical lattice; in this case, the final 
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FIG. 6: (color online) Profile of the condensate wavefunction 
at t = Is for the gap soliton generation shown in Fig. |3 
Dotted line shows the profile of the optical lattice. 



FIG. 5: (color online) Top: Temporal evolution (from left to 
right) of the axial density of the condensate in the case of 
adiabatic reduction of the lattice potential from e = S-Erccoii 
to e = -Erccoii during first t — 500ms. Middle: Evolution 
shown in Fourier space. Bottom: Temporal variation of the 
depth of the optical lattice. 



state corresponding to a gap soliton is composed by a few 
hundred of atoms. Moreover, the matter-wave gap soli- 
tons created in this way are very narrow, and they occupy 
only a few lattice sites. In addition, the loading efficiency 
is relatively low. For that reason, below we suggest and 
discuss in detail the other method for generating matter- 
wave gap solitons. In this second scenario, we propose 
replacing suddenly the longitudinal harmonic trap with 
a deep lattice, and then gradually reduce the lattice depth 
to a certain final value. This kind of the condensate dy- 
namics does not leave atoms as much freedom to escape 
the central region as they had before, ergo we can achieve 
higher values of the atomic densities in the soliton. More- 
over, by applying this method we can end up with very 
shallow traps, much below the depth accessible in the 
first scenario. In Fig. |S1 we present the result of our nu- 
merical simulations of the one-dimensional model for an 
adiabatic reduction of the lattice depth from e = SSrocoii 
to e = i^rocoii- During this dynamics, the size of the con- 
densate grows considerably due to a strong reduction of 
the lattice depth. In simulations shown in Fig. (SJ we were 
able to decrease the lattice depth to a level 2.5 times lower 
than in the first method discussed above, and we still ob- 
tained a stable state with about 500 atoms. The profile of 
the corresponding matter-wave gap soliton is presented 
in Fig. ini and the noticeable difference between Fig. |2| 
and Fig.iniis the width of the created gap soliton. Both 
gap solitons contain the same total number of atoms but. 



since in the second method we use much weaker lattice, 
the atoms spread over much larger region. 

To compare the efficiency of these two generation 
schemes, we perform additional numerical simulations 
and study how a variation of the initial number of atoms 
and the lattice depth influence the final number of atoms 
in the matter- wave gap soliton. In Fig.d we monitor the 
efficiency of the soliton generation by counting the atoms 
that remain after t = 0.5s within the window of 200 mi- 
crons around the central peak. This choice is somehow 
arbitrary, but the result does not depend qualitatively 
on the size of the window. We observe that our adia- 
batic scheme allows to achieve almost 100% generation 
efflciency in a broad range of the initial atom number. 
We also reveal that there exists a limit in the final num- 
ber of atoms that we can trap into the gap soliton. Both 
the schemes display saturation above 1000 atoms. We 
can estimate the maximal number of atoms that can be 
trapped in the soliton by employing the effective-mass ap- 
proximation and assuming that the size of the soliton can 
not be smaller than one lattice well. Using the parame- 
ters given above we estimate the upper limit to be around 
1200 atoms, that is by almost five times larger that that 
reported in the first experimental observation 12] . 

It is clear that the adiabatic scheme enables to achieve 
higher final population, and Fig. |S1 shows clearly the ad- 
vantage of this method. In this figure, we plot the final 
number of atoms in the gap soliton (we use the similar 
criterion as discussed above, starting with 500 atoms) as 
a function of the final value of the optical lattice depth. 
The adiabatic method leads to much smaller loss of atoms 
during the loading process, and it can be used to gener- 
ate matter-wave gap solitons in much shallower lattices. 
Solitons generated in this way span over quite large num- 
ber of lattice sites, and they exhibit much higher mobil- 
ity which can be useful for potential applications, see 
[l8l | . The latter feature that can be seen in Fig. |S1 which 
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Initial number of atoms 



FIG. 7: (color online) Number of atoms in the window 
[— 100/im, 100/im], which is approximately equal to the num- 
ber of atoms in the generated gap soliton, at t = 500ms, 
as a function of the initial number of atoms. Dashed line 
corresponds to the scenario with a constant lattice depth of 
e = 2.3-E'recoii, and the solid line- to the case of the adiabatic 
reduction of the lattice potential after switching it on from 
e = 8-Erecoii to £ = 2.3-Erecoii duriug the first 200ms. Straight 
line marks the available number of atoms. 
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FIG. 8: (color online) Number of atoms in the window 
[—100/im, 100/im] at t = 5s which is approximately equal to 
the number of atoms in the generated gap soliton as a func- 
tion of the final depth of the optical lattice. Initially, there 
are 500 atoms in the condensate. Dashed line corresponds to 
the case of a constant lattice depth, and the solid line- to the 
case of the adiabatic reduction of the lattice potential from 
£ = S-Erccoii to the final depth during the first 500 ms. 



shows a clear threshold for the soUton generation better 
pronounced for the adiabatic approach. 

Finally, we comment briefly on the soliton lifetime. In 
our model and numerical simulations we neglect inelastic 
losses, so that the timescale at which we follow the soliton 
generation dynamics remains quite long. Still, we have 
checked that including inelastic losses does not change 
significantly our main results, due to low densities of the 
localized structures generated, but it remains the main 
factor that will limit the lifetime of the gap solitons in 
real experiments. 



IV. COUPLED-MODE THEORY OF 
SELF-TRAPPING 



Fourier transform of the one-dimensional rescaled model 
'^^^ - y*(fc,i) - - 27r,t) -f ^{k + 27T,t)] + 



9_ 

'2t: 



/■ r + oc 

/ / ^{ki,t)-9*{k2,t)-^{k- ki+ k2,t)dkidk2. 



For the first Brillouin zone, i.e. \k'\ < tt, we define the 
Bloch functions u{k, k') of energy E{k') satisfying 

E{k')u(k,k') = -j'^ik, k'y-[u{k-2TT, k')+u{k+2TT, k')], 

with normalization condition J u{k,k[)u*{k,k'2)dk = 
(5(fc^ — ^2). According to our assumption mentioned above 
we can decompose the wavefunction 5'(fc,i) into these 
states 



In this section, we introduce a simple analytical model 
which allows to get a deeper insight into the physics of the 
soliton generation and self-trapping in the nonlinear sys- 
tems with repulsive interaction. Such a model describes 
the dynamics in the momentum space, but it is restricted 
by the first Brillouin zone of the bandgap spectrum. 

To derive the model, we decompose the condensate 
wavefunction in the basis of the Bloch functions and then 
show the spontaneous migration towards the region of the 
negative effective mass. 

Before the derivation, we introduce the dimensionless 
variables z' — z/d, t' — tuj, ^' — "^^Jd/N, where uj = 
h/{md'^) and N is the number of atoms, and define g = 
AndNaiu and e' — e/{huj). 

First, we omit the primes for simplicity and apply the 



^ik,t) = / <i>{k',t)u{k,k')dk', 



and obtain the equation for the amplitudes 
.a$(fc',t) 



i—^=E{k')^k\t) 



2tt 



-f 00 



<^>{k[ , t)$* {k'2,t)<i>{k'^,t)dk[dk'2dk'3 X 
u*{ko,k')u{ki,k[)u*{k2,k!2) x 



X u{kQ ^ ki + k2,k'2)dkQdkidk2. 

For a weak lattice we assume that u{k, k') « 5{k—k'); this 
approximation works well everywhere except the regions 
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close to the band edge. Using this assumption, we obtain 
.9$(fc,t) 



dt 



E{k)^k,t) 



(3) 



+ 



2tt 



We further simphfy the dynamics by introducing a two- 
mode approximation. Namely, we assume that to the 
first order our wavefunction consists of the constant and 
modulated parts, ^{k,t) = [a{t) + V2b{t) coa{k)] / ^/27r , 
E{k) = Eo + Eicos{k), where Ei < and jap + |6|2 = 1. 
Equation Q yields 



da 


i 


'dt ~ 


"71 


db 


i 


di ~ 


"71 



E,a--g\b\%. 



(4) 



After substituting 

ait) = 
bit) = 



1 - z(t) 



1 + zit) 



exp[i0o(t)]. 



exp[i(0o(i) - ^it))], 



we obtain the asymmetric Josephson junction equations 



d^ 

m 

dz 
'dt 



: COS( 



where we introduce A = 3(7/(4\/2|i?i |), and use rescaled 
time tnew = V2\Ei\t. To describe the migration of waves 
from the region of normal to anomalous diffraction, we in- 
troduce the quantity M(t), which stands for a difference 
between the population of atoms in the normal (central) 
and anomalous (outer) regions of the Brillouin zone 



Mit) 



it/2 /•-7r/2 

l$|2dfc-(/ m^dk- 

-7r/2 



/2 



If A/ > 0, the condensate is composed mostly of the 
matter waves from the normal diffraction region, and 
for M < - of the matter waves from the anomalous 
diffraction region. The temporal derivative of Al can be 
expressed as 



dM 

dt 



2V2 



A 



l\ dz 
'3) 'dt' 



(5) 



To apply this equation for describing the dynamics of 
the soliton generation, we assume that at f = the con- 
densate is in the ground state of a single lattice site. In 
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FIG. 9: (color online) Evolution of the parameter M(t) above 
the threshold for soliton formation (N = 500 atoms), cor- 
responding to Fig. 0^c) . Solid line shows the results of the 
one- dimensional GPE model, and dashed line shows results 
of the coupled-mode model, Eqs. (0J. 



M(t) 
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FIG. 10: (color online) Same as in Fig. El but in the case of 
N = 200 atoms, see Fig.IHb). 



this state, the first band is populated approximately uni- 
formly with a constant phase, so that a = 1 and b — 
and, consequently, z = —1 and Af (0) — 0. Since z cannot 
be smaller than — l,z will grow in time and from Eq. (j^J 
we see that M(t) will initially decrease for a repulsive 
condensate ( A > 0) . Hence the matter waves will change 
their state migrating to the regions of the negative effec- 
tive mass, allowing the gap soliton to be formed. 

This general picture has been confirmed by direct 
numerical simulations of both the full Gross-Pitaevskii 
equation and our simplified coupled-mode model, as 
shown in Fig. |51 and Fig. ^1 In Fig. |51 we present a 
comparison of the parameter Mit) calculated from one- 
dimensional GPE and Eq. Q for the case when the soli- 
ton is formed (500 atoms in the initial condensate) and 
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Fig.llOlcorresponds to the case below the threshold of the 
gap-soliton formation. In both cases the initial dynamics 
is similar and well described by our simple model. How- 
ever, the nonlinear dynamics observed in the simulations 
of the full model is more complicated because the energy 
exchange between the two major modes is accompanied 
by excitation of other modes in the system not taken 
into account by our simplified two-mode approximation. 
In the results obtained for the full model, this effect is 
observed as an effective damping of the oscillations. 

V. CONCLUSIONS 

We have revealed that the generation of matter-wave 
gap solitons can be easier than expected. We propose 
and demonstrate numerically two generation schemes in 
which a robust, long-lived stationary wavepacket in the 
form of a matter-wave gap soliton is created in a repulsive 
Bose-Einstein condensate placed into a one-dimensional 



optical lattice. The suggested generation method looks 
simple and efhcient, and it relies on a relaxation of the 
initial distribution of atoms to the appropriate soliton 
state. For this scheme, lifetime of the final state is limited 
only by the lifetime of the condensate. We presented a 
simple theoretical model to illustrate how the generation 
method works. 
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